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Multifrequency study of SNR J0533-7202, a new 
supernova remnant in the LMC 
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ABSTRACT 

We present a detailed study of Australia Telescope Compact Array (ATCA) ob- 
servations of a newly discovered Large Magellanic Cloud (LMC) supernova remnant 
(SNR), SNR J0533-7202. This object follows a horseshoe morphology, with a size 37 pc 
X 28 pc (1-pc uncertainty in each direction). It exhibits a radio spectrum with the 
intrinsic synchrotron spectral index of a=-0.47±0.06 between 73 and 6 cm. We report 
detections of regions showing moderately high fractional polarisation at 6 cm, with a 
peak value of 36±6% and a mean fractional polarisation of 12±7%. We also estimate 
an average rotation measure across the remnant of -591 rad m~'^. The current lack 
of deep X-ray observation precludes any conclusion about high-energy emission from 
the remnant. The association with an old stellar population favours a thermonuclear 
supernova origin of the remnant. 

Key words: supernova remnants - Large Magellanic Cloud - SNR 0533-7202. 



1 INTRODUCTION 

Supernova remnants (SNRs) are responsible for the distri- 
bution of heavy elements in the universe, influencing the 
chemical composition of the next generation of stars. SNRs 
have a signiflcant effect on their environment, heating up 
the surrounding gas and dust as the shock waves of the su- 
pernova explosion pass through. In turn however, SNRs are 
heavily impacted by their environment, as their evolution, 
structure and expansion is greatly affected by the density 
of the surrounding ISM. In the radio-continuum, SNR emis- 
sion is predominately non-thermal and will typically display 
a radio spectral index of a ~ -0.5 (defined by 5 oc u°'). 
However, this may significantly vary due to environmental 
factors, different stages of evolution and various types of 
SNRs. 

The Large Magellanic Cloud (LMC) is an irregular 
dwarf galaxy in close proximi ty to our own M ilky Way 
galaxy at a distance of 50 kpc (jMacri et al.l 120061 '). As a re- 



sult of this relatively close distance, objects in the LMC can 
be observed with significantly better spatial resolution and 
sensitivity than those in other galaxies. However, it is far 
away enough that we are able to assume all objects that 
lie within the galaxy are at very similar distances, therefore 
making estimates for various analysis methodologies such as 
extent and surface brightness more accurate, as these val- 
ues depend on accurate distances to the object. In contrast, 
objects in the Milky Way can be hard to find accurate dis- 
tances for, and therefore leads to imprecise measurements. 
The LMC is also a desirable environment for various astro- 
nomical studies as it contains some of the most active star 
forming regions in our Local Group of galaxies as well as 
residing outside of th e Galactic plane, at a mode rate in- 
clination angle of 35° (|van der Marel fc Cionill200ll ). which 
greatly minimises the interference from stars, gas, and dust 
within the LMC. 

In this paper, we report on a newly discovered LMC 
SNR, SNR J0533-7202. Radio observations of this object 
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were also taken by iFilipovic et alj (|l995l . Il998al fbr) in their 
surveys of the Magellanic Clouds, however, they did not 
classify this object. The new observations, data reduction 
and imaging techniques are described in Section 2. The 
astrophysical interpretation of newly obtained moderate- 
resolution total intensity and polarimetric images in com- 
bination with the existing Magellanic Cloud Emission Line 
Survey (MCELS) images are discussed in Section 3. 



2 OBSERVATIONS AND DATA REDUCTION 

We observed SNR J0533-7202 with the ATCA on the 15*'' 
and IG**" of November 2011, using the new Compact Array 
Broadband Backend (CABB) at array configuration EW367 
and at wavelengths of 3 and 6 cm (i/^9000 and 5500 MHz). 
Baselines formed with the 6*'' ATCA antenna were excluded, 
as the other five antennas were arranged in a compact con- 
figuration. The observations were carried out in the so called 
"snap-shot" mode, totaling ~50 minutes of integration over 
a 14 hour period. PKS B 1934-638 was used for flux den- 
sity calibratioio and PKS B0530-727 was used for secondary 
(phase) calibration. The phase calibrator was observed twice 
every hour for a total 78 min utes over the whole observ- 
ing session. T he miriaeo (Sa ult et al.l 119951 ) and karma 
ICoochl I1995I ) software packages were used for reduction 
and analysis. More information on the observing procedure 
and other sour c es obs erved i n this project can be found in 
iBozzetto et"ai] (|2012al lbh and lde Hortaet al.l l|2012l '). 

The CABB 2 GHz bandwidth is a 16 times improve- 
ment from the previous 128 MHz, and with the new higher 
data sampling has increased the sensitivity of the ATCA by 
a factor of 4. The 2 GHz bandwidth not only aids in high 
sensitivity observations, but also allows data to be split into 
channels which can then be used for measuring Faraday rota- 
tion across the entire bandwidth, at frequencies close enough 
that the n x 180° ambiguities prevalent when making an es- 
timate between distant frequencies, are no longer an issue. 

I mages were f ormed using miriad multi-frequency syn- 
thesis ISault fc Wieringa 1994 ) and natural weighting. They 
were deconvolved with primary beam correction applied. 
The same procedure was used for both U and Q stokes 
parameter maps. 

The 3 cm image (Fig.[T]) has a resolution (full width half 
maximum (FWHM)) of 21.6" x 15.0" (PA=47.2°). Similarly, 
we made an image of SNR J0533-7202 at 6 cm (seen as 
contours in Fig. [TJ which has a FWHM of 34.1"x26.1" at 
PA=45.5° and an estimated r.m.s. noise of 0.3 mjy/beam. 



3 RESULTS AND DISCUSSION 

This new LMC remnant exhibits a horseshoe mor- 
phology (Fig. [1]), centered at RA(J2000)=5''33'"46.5", 
DEC(J2000)=-72°02'59". We selected a one-dimensional in- 
tensity profile across the approximate major (NW-SE) and 
minor (NE-SW) axis (PA=45°) (Fig.[2l) at the 3cr noise level 
(0.9 mjy) to estimate the spatial extent of SNR J0533-7202. 



Its size at 6 cm is 152" x 115" with a 4" uncertainty in 
each direction (37x28 pc with a 1 pc uncertainty in each 
direction). We did not detect any [Oni] or [Sll] emission 
in the Magellanic Cloud Emission Line Survey (MCELS) 
(jSmith et al.l 120061 ) images. However, there is some tenta- 
tive and very faint Hq emission possibly associated with 
this SNR which wasn't evident anywhere else in the field 
surrounding the remnant. Although, more sensitive obser- 
vations are needed. 

An X-ray source at the rim of the SNR was detected 
in the course of the ROSAT all-sky su rvey and was give n 
the identifier IRXS J053353. 6-720404 (|Voges et al.lll999l ). 
However, the likelihood of existenc4J in these observations 
was only 7, meaning only a 3.3cr detection, and therefore 
not much could be done with respect to the X-ray analy- 
sis, apart from plotting the ROSAT position (Fig. [T]). This 
ROSAT position is slightly south-west of one of the radio- 
bright regions of the remnant. Low statistics preclude any 
classification between an extended or compact source. 

This S NR did not appear in the Spitzer mosaic images 
of the LMC JMeixner et al.ll2006l ). neither at the 3.6, 4.5, 5.8, 
and 8 /xm wavelengths of the IRAC instrument (JFazio et al.l 
l2004h . nor i n the 24, 70, and 160 /^m bands of the MIPS 
instrument (jRieke et al.ll2004l ). suggesting that there is no 
association with mid or near-infrared wavelengths. There are 
no OB star candidates in the Magellanic Clouds Ph otometric 
Survey (MCPS) catalogue jZaritskv et al.l |2004| ) within a 
100 pc radius around the centre of the remnant, and the 
star formation history map of the LMC l|Harris fc Zaritskvl 
12009) shows no recent episode of star formation activity in 
the neighbourhood. The association of SNR J0533-7202 with 
an old stellar population favours a thermonuclear supernova 
origin of the remnant. 

We based the spectral energy distribution (SED) on our 
own integrat ed flux estim a tes, c oupled with the 73 cm mea- 
surement bv lClarke et al.l (|l97y). These values are shown in 
Table 1 and then used to produce a spectral index graph 
(Fig. El). The point source (ATCA J0534-7201; see Fig. [TJ 
to the east of SNR J0533-7202 was unresolved from the 
SNR in the 73 cm survey (MC4; IClarke et all 1 19761 ). We 
estimate 36 cm flux density measurements from the Molon- 
glo Synthesis Telescope (MOST) mosaic image (as described 
i n Mills et al 1984 i) and a 36 cm SUMMS mosaic image 
(jMauch et al.ll2008l ). We point that these values differ some 
25%, which is most likely due to the higher sensitivity of 
the MOST image and its greater UV coverage. The 20 cm 
flu x density was measured from a mosaic image published 
bv lHughes et al.l (|2007i ). Two different sets of images were 
used to estimate integrated flux densities at wavelengths of 
6 & 3 cm. The first set (5500 & 9000 MHz) contains our 
CAB B observations mer ged with the mosaic visibility files 
from iDickel et al.l (|2005h . These observations use the EW 
367 array, which has a shortest baseline of 46 m and as a 
result, has missing flux from the lack of short spacings. The 
second set of images (4800 fc 8640 MHz) are from the mo- 
saic images published by Dickel et al. (2010), which used 
EW 367 and EW 352 arrays at both frequencies as well as 
Parkes at 4800 MHz but not at 8640 MHz, as the Parkes 



^ Flux densities were assumed to be 5.098 Jy at 6 cm and 2.736 
at 3 cm. 

^ http://www.atnf.csiro.au/computing/softwarc/miriad/ 



•^ In the ROSAT source detection a likelihood L is associated to a 
probability of a detected source being real P by P = 1 - exp(-L). 
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survey at 3 cm did not extend that far south. This survey 
data had more spatial frequency coverage, particularly for 
the extended emission, but less sensitivity. We can see the 
detrimental effect of missing short spacings in the 3 cm flux 
density measurements, where they fall well below the trend 
of the SED at higher wavelengths (Fig.[31). Due to the signif- 
icant impact of the missing shorter spacings (and as a result, 
missing flux) , we omit the 3 cm measurements from the cal- 
culation, leaving all the frequencies up to 6 cm (5500 MHz), 
and thus a spectral index of q:=-0.47±0.06. 

Fractional polarisation (P) was calculated at 6 cm us- 
ing: 



Si 



(1) 



where Sq, Su and Si are integrated intensities for Q, U and 
/ Stokes parameters (Fig. [4]). Our estimated peak value is 
36±6% with a mean fractional polarisation of 12±7%. This 
level of fractional polarisation is relativel y high when com- 
pared to various othe r SNRs in the LMC JCaiko et al.ll2009l : 
ICrawford et all |201G| : iBozzetto et all l2012al ') and would be 
(theoretically) expected for an SNR with a radio spectrum 
of around or less than -0.5 (Rolfs & Wilson 2003). This may 
indicate varied dynamics along the shell. 

Polarisation position angles were taken from across the 
2 GHz bandwidth (at 5500 MHz split into 128 MHz chan- 
nels) and used to estimate the Faraday rotation for this 
SNR. The result of this can be seen in Fig. [S] with the 
open boxes representing negative values of rotation measure 
and the filled in boxes representing positive rotation mea- 
sure. The rotation measure varies quite significantly along 
the SNR with the most change being negative rotation mea- 
sure near the peak intensity in the western region of the 
remnant. The average rotation measure across the SNR was 
-591 rad m~^. This value is n early double that of the plerion 
in SNR G326.3-1.8 (JDickel e t al. 2000) and significantly ex- 
ceeds what is typical of ' large' values in the LM C and Milky 
Way of ±250 rad m"^ (|Dickel fc Milnelll995h . It would be 
best to treat this value with some level of caution, as at 
higher radio frequencies - such as our 6 cm observations 
- the amount of rotation measure expected from a SNR is 
within the same range as the expected error. To mitigate this 
error and achieve a more reliable value of rotation measure, 
additional observations would need to be taken, preferably 
at lower radio frequencies where we would expect a higher 
level of rotation measure. 

We were also able to estimate the magnetic field 
strength f or this SNR based on the equipartition formula 
Arbutina et al.l l|2012f ). This formula is based 



the lBelll (119781 ) diffuse shock acceleration (DSA) theory. By 
using the spectral index value q=-0.5, we get an equiparti- 
tion value for SNR J0533-7202 of ~45 jxC with an estimated 
minimum energy of Emin ~ 9.4x10*® ergs. 

From the position of SNR J0533-7202 at the surface 
brightness to diameter (E - D) diagram (( D, S) = (32.5 pc, 
6.2 X 10"^^ W m-2 Hz-i sr'^)) by iBerezhko fc Volkl 
l|2004h , we can estimate that SNR J0533-7202 is likely to 
be an SNR in the late energy conserving phase, with an ex- 
plosion energy between 0.25 and 1 x 10^^ ergs, which evolves 
in an environment of density ~1 cm~'^. 
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Figure 1. ATCA observations of SNR J0533-7202 at 3 cm 
(9000 MHz) overlaid with 6 cm (5500 MHz) contours. The con- 
tours are 3, 6, 9, 12 & I5cr. The ellipse in the lower left corner 
represents the synthesised beamwidth (at 6 cm) of 34.1 "x26.1". 
The sidebar quantifies the pixel map. The overlaid circle shows 
the position of a weak X-ray source seen by ROSAT. 



4 CONCLUSION 

We have added a new SNR to the LMC SNR population 
through conducting a high resolution radio-continuum 
study of SNR J0533-7202. We report a relatively large 
SNR with an extent of ~152"xll5" (~37x28 pc), and a 
radio spectral index with q=-0.47 between 73 and 6 cm. 
We estimate fractional polarisation of the remnant at 6 cm 
with a peak of 36±6% and a mean integrated value of 
12±7%. The lack of recent star formation activity around 
the remnant makes a thermonuclear supernova origin more 
likely. 
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Figure 2. The top image sliows the 6 cm intensity image of 
overlaid with the approximate major (NW-SE) and minor (NE- 
SW) axis. The middle and lower images show the 1-dimensional 
cross-section along the overlaid lines in the top image, with a 
superimposed line at Scr. 
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Figure 3. Radio-continuum spectrum of SNR J0533— 7202 and 
the point source ATCA J0534-7201. The pink dotted line rep- 
resents the spectral index of the SNR and the blue dashed line 
represents the spectrum of the point source. 
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Figure 4. B-field vectors overlaid on 6 cm contours (3, 6, 9, 12 
and 15(t) of SNR J0533-7202 from ATCA observations. The el- 
lipse in the lower left corner represents the synthesised beamwidth 
of 34.3" X 26.0" and the line below the ellipse shows a polarisation 
vector of 100%. 
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Figure 5. Faraday rotation measure of SNR J0533-7202 over- 
laid on 6 cm (128 MHz bandwidtii) ATCA contours (3, 8 & 
13o"). Filled squares represent positive rotation measure while 
open squares represent negative rotation measure. The eUipse 
in the lower left corner represents the synthesised beamwidth of 
46.7" X 27.9" and the box in the lower right represents a rota- 
tion measure of 1000 rad/m-^. The width of the boxes scale with 
rotation measure. 
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